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ABSTRACT 

We present new dynamical models of the merger remnant NGC 7252 which include star 
formation simulated according to various phenomenological rules. By using interactive 
software to match our model with the observed morphology and gas velocity field, we 
obtain a consistent dynamical model for NGC 7252. In our models, this proto-elliptical 
galaxy formed by the merger of two similar gas-rich disk galaxies which fell together 
with an initial pericentric separation of ~ 2 disk scale lengths approximately 620 Myr 
ago. Results from two different star formation rules — density-dependent and shock- 
induced — show significant differences in star formation during and after the first 
passage. Shock- induced star formation yields a prompt and wide-spread starburst at 
the time of first passage, while density-dependent star formation predicts a more slowly 
rising and centrally concentrated starburst. A comparison of the distributions and ages 
of observed clusters with results of our simulations favors shock-induced mechanism of 
star formation in NGC 7252. We also present simulated color images of our model of 
NGC 7252, constructed by incorporating population synthesis with radiative transfer 
and dust attenuation. Overall the predicted magnitudes and colors of the models are 
consistent with observations, although the simulated tails are fainter and redder than 
observed. We suggest that a lack of star formation in the tails, reflected by the redder 
colors, is due to an incomplete description of star formation in our models rather than 
insufficient gas in the tails. 

Key words: galaxies: formation - galaxies: individual (NGC 7252) - galaxies: inter- 
actions - galaxies: kinematics and dynamics - galaxies: structure 



1 INTRODUCTION 



NGC 7252 is a well-studied merger remnant (|Toomrelll977l ; 
ISchweizer! [l983h . which provides strong observational evi- 
dence supporting the idea that at least some disk galaxy 
mergers produce elliptical galaxies. Its luminosity profile 



follows a 



rV4 



law (ISchweizer 19821), as do the profiles o f 



elliptical galaxies (|de Vaucouleursl Il953l ; iKormandvl 1 19771 ) . 
The core contains a central counter-rotating disk of ion- 
ized gas within 8" 2.5 kpcQ) of the center, providing 
a kin ematic signature of a major merger event in the recent 
past (!Schweizerlll982T ). Molecular gas is also closely associ- 



* E-mail: chien@stsci.edu (LHC); barnes@ifa.hawaii.edu (JEB) 
1 In order to compare with previous observational results, 
throughout this paper we adopt Ho = lOO/i = 75 km s _1 Mpc -1 . 
This places NGC 7252 at a distance of 64.4 Mpc and a projected 
scale of 1" = 312 pc. 



ated with the central disk (|Dupraz et al.lll990l : IWang et al.l 
1992). The two tidal tails are rich in HI, s uggestin g that 
the progenitors were late-type spirals ([Hibbard et al.lll994l 
hereafter HGvGS94). About 500 candidate young clusters 
are foun d in the main body in the H ubble Space Telescope 
images (|Whitmore et al.lll993l : iMiller et alJll997h . indicat- 
ing that clusters can form in large numbers in the violent 
star formation that accompanies a galaxy merger, and that 
cluster formation can increase the specific globular cluster 
frequency of the merger remnant. Moreover, since the clus- 
ters are single burst stellar populations, the gravitationally 
bound clusters surviving several 10 8 years provide a record 
of the merger's star formation history. Ages and metallici- 
ties of these clusters can be used to reconstruct that history; 
ISchweizer fe Seitzed (|l998l . hereafter SS98) found five young 
globular clusters with ages of 400-600 Myr, which presum- 
ably formed shortly after the first close encounter of the 
two gas-rich disk galaxies that later merged to produce the 
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present-day remnant. Metallicities of these clusters support 
the idea that elliptical galaxies with bimodal globular cluster 
systems can form through major mergers. 

In addition to the extensive observations of NGC 7252, 
many numerical models have been made to inve stigate the 
history of the system. iHibbard fe Mihod (1995, hereafter 
HM9 5) revised encounter pa r ameters from prev ious stud- 
ies (jBorne fe Richstond Il99ll : iMihos et al.1 1 19931 ) with de- 
tailed Hi observations from the Very Large Array (VLA, 
HGvGS94). Their dynamical simulation determined the age 
of NGC 7252 to be about 773 (~ 580/T 1 ) Myr after peri- 
apse of the initial orbit an d predicted the fut ure infall rate 
of the tidal tail materials. Mi hos et al.l (|l998h further con- 
strained the dark halos in the NGC 7252 progenitor galax- 
ies. These models matched the morphology and kinematics 
of the system very well, and have provided robust results 
for the dynamical history. However, gas dynamics and star 
formation have not yet been included in any dynamical mod - 
eling of NGC 7252. lFritze-v. Alvensleben fe Gerhardl (|l994h 
attempted to reconstruct NGC 7252 's star formation history 
and spectro-photometric properties using an evolutionary 
spectral synthesis model, but were forced to postulate that 
the interaction-induced star formation began well before the 
first passage. 

In this paper we use simulations including gas dynamics 
and star formation to construct possibl e star formation his- 
tories for NGC 7252. We follow iBarnesI ((2004) in examining 
how two different star formation rules — density- dependent 
and shock-induced — affect the star formation history. In 
addition to matching the observed morphology and H I kine- 
matics, we compare our star formation models with the ages 
of young globular clusters observed by SS98. We examine the 
different star formation models and try to constrain the trig- 
gers of star formation in merging galaxies. Moreover, we also 
compare the modeled colors of NGC 7252 with the observed 
values. Our goal is to match the overall morphology, kine- 
matics and photometric properties while insuring that the 
initial disks have reasonable parameters for late-type spi- 
rals. This represents a new and significant step beyond pure 
dynamical modeling. This paper is organized as follows: in 
Section 2 we describe the methodology and the simulation 
technique, while in Section 3 we describe the matching pro- 
cess to the HI observation as well as our best-fit model. In 
Section 4 we present the predicted star formation history 
and a comparison with the ages of SS98 clusters. We discuss 
the simulated images of NGC 7252 and a comparison of the 
photometric properties in Section 5, and finally give conclu- 
sions in Section 6. For the purpose of a consistent discussion 
throughout the paper, we introduce physical scales in Sec- 
tion 2 using the scaling factors of the best-fit model, which 
will be described in Section 3. 



2 NUMERICAL TECHNIQUES 
2.1 Galaxy Models 

We use bul ge /disk /halo galaxy mode ls similar to those de- 
veloped bv Barnes fc Hibbardl (|2009l . Ap pendix B) Eac h 
galaxy has a bulge with a shallow cusp (|Hernauistl[l99Qh . 
an exponential disk with constant scale height, and a trun- 
cate d dark matter halo with a Navarro- Frenk- White pro- 
file ([Navarro et al.l Il997l ) . The mass and scale lengths of 



Table 1. Parameters of the components of each progenitor 
galaxy. 



Component 


Value 




0.62xl0 10 M 


M-disfc, stellar 


1.40xl0 10 M 


~^-disk, gas 


0.47X10 10 M 


M hato 


9.96xl0 10 M 


^bulge 


0.44 kpc 


^disk, stellar 


0.28 kpc 


z disk, gas 


0.11 kpc 


T> C 

r^disk 


1.84 kpc 


&halo 


5.53 kpc 



a Bulge scale length. 
b Disk scale height. 
c Disk radial scale length. 
d Halo scale radius. 



each comp onent is listed in Table [TJ The stellar disk's Q 
parameter (jToomrd 1 19641 ) is initially greater than 1.45 at 
all radii, insuring local stability and helping to suppress 
spontaneous bar formation. The total mass of the disk is 
Mdiak - 1-87 x 10 10 M , slightly smaller than the Milky 
Way but reasonable for typical spiral galaxies. The disk 
has a circular velocity v c (R = 3Rdisk) — 176 km s _1 . We 
also estimate its maxi mum rotation velocity , v m ax , from the 
Tully-Fisher relation (|Tullv &; Fisher! 1 19771 ). assuming that 
the total mass, M, is proportional to the fourth power of the 
circular velocity, or M oc . Using the Milky Way as a refer- 
ence (i.e. a disk mass of Mmw — 6.0 x 10 10 Mq and circular 
velocity of 220 km s _1 ), we predict a circular velocity for our 
galaxy model v max (M disk /M M w) 1/4 x 220 164 krns" 1 . 
This value is slightly smaller than the circular velocity of 
our model galaxy, but well within the observed scatter of 
the Tully-Fisher relation. 

The total halo to disk-plus-bulge mass ratio of our 
galaxy is 4:1. Thus our galaxy models have relatively low- 
mass halos compared to those expe cted in som e scenarios 
of galaxy formation in CDM (e.g. IMq et~ai1ll998h . although 
they are not very different from the models used in previous 
simulations of NGC 7252. We conjecture that the orbital de- 
cay of a pair of galaxies following a relatively close passage 
is not very sensitive to the outer halo structure; this could 
be explored in experiments embedding galaxy encounters in 
a self-consistent cosmological context, but it is beyond the 
scope of the present study. 

Gas is initially distributed like the stellar disk and 
follows an isothermal equation of state with temperature 
T ~ 10 4 K. The star formation history depends on the 
amount of gas that is available in both galaxies before the 
encounter. SS98 suggested that NGC 7252 is formed by 
the merger of two gas-rich disk galaxies, possibly of Hub- 
ble type Sc. Accordingly, we made the progenitor galaxies 
gas rich, with gas comprising 25% of the disk mass. The 
total number of particles in each galaxy is 108544, broken 
down into 8192 bulge particles, 18432 disk particles, 32768 
halo particles, and 49152 gas particles, which are evolved us- 
ing a smoothed particle hydrodynamics (SPH) code. Grav- 
itational forces calculated by a modified N-body tree code 
with a Plummer smoothing of e ~ 0.17 kpc. 
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2.2 Star Formation Prescriptions 

The star formation rules incorporated in our simulation fol- 
low iBarnesI ([2004) . We adopt the star formation prescrip- 
tion: 

* Pg U , (1) 

where the star formation rate, p*, depends on both gas den- 
sity, p g , and local rate of mechanical heating due to shocks, 
ii. T his prescripti on is implemented in a probabilistic fash- 
ion (jBarnesI 12004). At each time step, the probability of a 
gas particle undergoing star formation is calculated accord- 
ing to the assumed star formation rule. If a gas particle 
satisfies the probabilistic criteria for star formation, it is 
instantly and completely transformed into a stellar parti- 
cle. This scheme provides a Monte- Carlo description of star 
formation- consistent with the overall Monte- Carlo nature 
of N-body simulation- while keeping the total number of 
particles fixed and avoiding complications associated with 
composite star/gas particles. For later analysis, each stellar 
particle is tagged with the value of time at its birth. 

We considered two limiting regimes for NGC 7252: 
density- dependent and shock-induced star formation. The 
former star formation model is solely depending on the ac- 
cumulation of gas, with n = 1.5 and m = mimicking a 
Schmidt law (jSchmidtl fl959l : iKennicuttl [l998h . The latter 
model uses n = 1.0 and m = 0.5, favoring star fo rmation 
i n reg ions corresponding to shocks. As described in IBarnesI 



(2004), the constant C* has to be adjusted to obtain equiva- 
lent results in simulations with different spatial resolutions, 
and to insure that different simulations consume about the 
same amount of gas. We adopt C* = 0.018 for the density- 
dependent model and C* = 0.5 for the shock-induced model. 

2.3 Star Formation Prehistory 

In the previous section we describe how we model ongoing 
star formation, but in order to calculate the photometric 
properties, we also need to specify the pre-encounter star 
formation history of the initial galaxies. We setup the pre- 
encounter star formation history of the initial galaxies as 
shown in Figure [1] First, we assume the bulge formed be- 
tween 11 and 10 Gyr before the start of the simulation (t 
= 0) ; since the integrated photometric properties of such an 
old population are insensitive to its detailed star formation 
history, we adopt a constant star formation rate (SFR) of 
6.2 Mq yr _1 during bulge formation. Second, we assume 
the disk began forming 10 Gyr before the start of the sim- 
ulation and adopt an exponentia l form for the disk's SFR 
(jSearle et al.lll973l : ISandagdfl986h . 



ty(t) = M e 



(2) 



where Mo (— 1.8 Mq yr _1 ) matches the dynamically- driven 
SFR at the start of the simulation. To determine the e- 
folding factor, a, we integrate equation J2} from t — —10 Gyr 
to t = and require the result to be equal to the mass of 
the stellar disk. In our simulation, each of the progenitor 
galaxies has a stellar disk mass of 1.40 x 10 10 M , so this 
factor is found to be a ~ —0.054 Gyr -1 . The SFR thus rises 
slowly throughout the previous 10 Gyr, with an average of 
1.4 M yr" 1 . This gives the characterization parameter b 
(Scalo 1986), defined as the ratio of the current SFR to the 




t (Gyr) 

Figure 1. Star formation prehistory of the progenitor galaxies 
(in Mq yr _1 ) as a function of time (in Gyr); t = is when the 
simulation starts. Disk populations began forming 10 Gyr ago and 
bulge populations formed 1 Gyr before the disk, corresponding to 
the spike. The small fluctuations in the star formation rate reflects 
our monte-carlo implementation of star formation. 



past SFR averaged over the age of the disk, a value of about 
1.3, consistent with typical late- type spirals. We emphasize 
that this star formation history is schematic; it is broadly 
consistent with known constraints, but other scenarios are 
possible. For example, instead of using a slightly unconven- 
tional e- folding factor a < 0, we could have given the disk 
a constant SFR of 1.8 Mq yr _1 and a maximum age of 
7.7 Gyr; this would make the the pre-encounter disks a bit 
more luminous but would have very little effect on the over- 
all outcome of our modeling. 

Finally, the birth-dates of individual bulge and disk par- 
ticles are assigned using a Monte-Carlo procedure. Bulge 
particle birth-dates are drawn from a uniform distribution 
between —11 and —10 Gyr, reflecting the constant SFR of 
the bulge assumed above. Disk particle birth-dates are as- 
signed by sampling the disk SFR \£(£) defined by equation 
(2). 



2.4 Modeling Photometric Properties 

The photometric properties of the modeled galaxies are de- 
termined by convolving the simulation output with a simple 
stellar population model. Using the conventional notation 
for SPH, the luminosity density, j, in waveband X, can be 
expressed as: 



jx(r,t) = ^2w(Y-Y i ,h)-tx(t-t i ), 



(3) 



where W is the kernel function, h is the smoothing length, 
and tx(r) is the luminosity of a stellar particle with age r. 
The sum runs over all stellar particles; particle i has position 
Yi and birthdate U. 

To evaluate part icle luminosities, we use 
iBruzual fe Charlotl (|2003l . hereafter BC03) m odels for 
a stel lar population with solar metallicity and IChabrierl 
(2003) initial mass function (IMF). One somewhat subtle 
issue concerns the mass returned to the interstellar medium 
(ISM) via winds, planetary nebulae, and supernova ejecta. 
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Following the normal population synthesis methodology, 
BC03 tabulate luminosities as functions of time for a stellar 
population formed in a burst of a fixed initial mass. The 
actual mass still in stars declines as matter is returned to 
the ISM; with the specific IMF and metallicity we have 
chosen, a stellar population loses ~ 12, ~ 28, ~ 39, and 
~ 49 percent of its mass after 10 Myr, 100 Myr, 1 Gyr, and 
10 Gyr, respectively. However, this process has not been 
included in the N-body simulation^]; particles representing 
stellar components have fixed masses, and the total number 
of such particles is either fixed (without star formation) 
or monotonically increasing (with star formation). The 
question is how to reconcile population synthesis, which 
includes mass loss, with N-body simulations, which do 
not. We chose to take the masses of our stellar particles 
"at face value" when assigning luminosities; we therefore 
correct BC03's tabulated luminosities upward so that 
they represent a population with a constant stellar mass. 
As indicated by the mass-loss figures given above, this 
correction is small for relatively young stellar particles 
but becomes significant for older ones. For example, when 
evaluating the luminosity of a bulge particle with an age 
of 10 Gyr, we multiply the tabulated luminosity by a 
factor of (1 - 0.49)" 1 ~ 1.96. Since each stellar particle, 
including those formed during the simulation, has a definite 
birth-date, it is straightforward to apply this correction on 
a particle- by-particle basis when assigning luminosities. 

The surface brightness along each line of sight, Ix, is 
then given by the standard equation of radiative transfer 
with no scattering: 



dlx 

dz 



: jx — pgKxIx, 



(4) 



where z is the distance traveled in the direction of the ray to 
the observer, p g is the density of the gas and kx is the opac- 
ity. We define the opacity according to the values of Milky 
Way so that a gas surface density of 2.0 x 10 7 Mq kpc -2 
yields an extinction of Ay — 1.00, A\j — 1.53, and Ab — 
1.32. Total magnitudes and colors are obtained by summing 
the surface brightness of the entire galaxy at each waveband. 



3 MORPHOLOGY AND KINEMATICS 

The orbital geometry of our simulations is largely based on 
HM95. The initial orbits were parabolic, with an pericen- 
tric separation of 2ARdi S k, or r p = 4.4 kpc. We found that 
the angles specifying the orientation of the progenitor disks 
in HM95 were slightly ambiguous; HM95's definition of the 
argu ment of pericente r , uj, d iffers from the definition given 
bv iToomre fe Toomrd (|l972h . Fortunately we were able to 
match the geometry of the initial conditions based on the 
proj ections on th e x-y and y-z planes (cf. Figure AVIII-1 
m lHibbardl fl995). obtaining inclination, z, and arguments, 
cj, of (inw ,ojnw) — (—40,90) and (ie,we) = (70,130), for 
the disks giving rise to the Northwest tail and the Eastern 
tail respectively. In the rest of this paper, this combination 
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2 Preliminary experiments by JEB suggest that a Monte-Carlo 
scheme for returning stellar particles to the gas phase is workable, 
but further development of this scheme is beyond the scope of this 
paper. 



Figure 2. A snapshot image of Identikit 1.0 ([Barnes fc Hibbardl 
2009). Top Left: Optical image of NGC 7252 overlaid with HI 
contours (HGvGS94). Top Right: Y — V z posit ion- velocity dia- 
gram. Bottom Left: V z — X position- velocity diagram (HM95). 



of pericentric separation and orbital angles is referred to as 
the HM95 setup for convenience. We also experimented with 
several combinations of different angles and pericentric sep- 
arati ons close to this se t up. 

iBarnes &; Hibba"rdl (|20Q9h have developed a software 
package, Identikit 1.0, which allows rapid exploration of the 
encounter parameter space to reproduce the observed mor- 
phology and kinematics of interacting disk galaxies. In the 
following we describe the steps of our matching process, 
which utilize this interactive software extensively. Figure [2] 
shows an image from the software, with the observed im- 
age and HI velocity data of NGC 7252 projected on each 
plane. The software color codes the particles according to 
their galaxy of origin, and displays them according to their 
projected velocities and positions. While the optical image 
and gas distribution show the observed morphology to be 
matched, the aligned panels of gas velocity provide strong 
constraints on the dynamics of the galaxy. 

First, an initial merger model was constructed with the 
HM95 setup, and self- consistent simulations incorporating 
both star formation rules were run based on this setup. For 
this simulation, we used the interactive display of Identikit 
to find a best match to the observed morphology and kine- 
matics. We then calculated the photometric properties of 
this model and compared them with the observed values, 
however we found that this model cannot explain the ages 
of young globular clusters observed by SS98, and neither 
the predicted star formation histories nor the photometric 
properties from this setup appeared consistent with the ob- 
servations. 

Thus the next step was to try different angles around 
the HM95 setu p as well as different peri centric separations. 
As described in Barnes fe Hibbardl (|2009L cf. their Figure 3), 
in addition to allowing the exploration of the parameter sub- 
space with a fixed set of initial conditions, Identikit also 
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Figure 3. The final best-fit model of NGC 7252, displayed as 
color-coded particles on top of the observed data, according to 
their galaxy of origin. Red particles relate to the progenitor galaxy 
that formed the E tail, while blue particles relate to the galaxy 
that formed the NW tail. Bottom Right: Modeled NGC 7252 
viewed from the projected X — Z plane, where Z is the direction 
along the line of sight. 



provides instantaneous access to the outcome of different 
encounter geometries. We used this capability by changing 
the pericentric separation, r p , examined Identikit simula- 
tions with different choices of (i, uj) around the HM95 setup, 
and found possible combinations of (i, uj) and viewing di- 
rections which matched the morphology and kinematics of 
NGC 7252. We then ran self-consistent simulations incorpo- 
rating both star formation rules based on this setup, and 
compared the star formation histories as well as the photo- 
metric properties to the observed values. If the results were 
inconsistent with the observations, this process was repeated 
with another choice of r p and the corresponding combina- 
tion of (z, uj), until a plausible match to the observational 
properties was found. 

The range of r p that we explored was 1.4 to 3Rdisk, 
or r p = 2.6 kpc to 5.5 kpc. The inclination angle inw was 
fixed as %nw — —40 in all the experiments, while the other 
angles were varied in the range of ujnw = 66 to 90, %e — 60 
to 70, and uje — 130 to 137. We emphasize here that, within 
the restrictions of time and machine power as well as the 
limitations of our numerical technique, an exact match to 
NGC 7252 is not possible; our goal is to obtain a star for- 
mation model that approximates the observational proper- 
ties and provides reasonable predictions of NGC 7252 's star 
formation history. In the following we present our final best- 
fit model and discuss some of the discrepancies between our 
model and the observations. 

The best-fit model assigns the two progenitor galaxies 
orbital angles of (inw^nw) = (—40,66) and (ie^oje) — 
(60, 137) and adopts a pericentric separation of r p = 
1.92 Rdisk — 3.5 kpc. The best match to the present occurs 
at c± 620 Myr after periapse and is shown in Figure [3] The 
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T (Myr) 

Figure 4. Global star formation history (in Mq yr _1 ) of 
NGC 7252 using two star formation rules. Periapse occurred at 
T ~ —620 Myr, indicated by the dashed line, and the best-fit 
time is at T = 0. Black and grey lines represent the SFR using 
the shock-induced and the density- dependent star formation rules 
individually. Horizontal lines at the bottom indicate the duration 
of each burst described in the text. 



overall morphology and kinematics of our simulation agree 
very well with the optical image and Hi data. In our best 
fit, the system is observed ~ 60° above the orbital plane 
and ~ 70° behind the line connecting the galaxies at orbital 
periapse. 

As in other simulations of gas-rich merger remnants 
(e.g. Barnes 2002), our model contains a rotating gas disk 
superficially sim ilar to the one observed in NGC 7252 
( Sc hweizerl Il982h . However the disk in our model appears 
nearly edge-on at the time of best match, of which ob- 
servations indicate an inclination of ~ 40°. The simulated 
disk seems to be precessing in the gravitational potential of 
the remnant, with a rate which presumably depends on the 
details of the model, including the initial mass profiles of 
the progenitor s and the smoothing of the force calculations 
(Barnes 2002). The orientation of the disk at the instant 
of best match might be improved by adjusting these details 
while preserving the overall choice of initial orbit, disk ori- 
entations and viewing angles; however this is beyond our 
resources. 

4 STAR FORMATION HISTORY OF NGC 7252 

4.1 Interaction-Induced Star Formation 

Figure H] shows the overall star formation histories predicted 
using the two different star formation models. In this section 
we denote T = as the best-fit time, i.e. the present; thus 
the periapse occurred at T c± —620 Myr. In the early stage of 
interaction, the two galaxies fall towards each other and in- 
terpenetrate at periapse. This produces a prompt starburst 
in the shock-induced simulation, where star formation quick- 
ens slightly before periapse, peaking at the time of closest 
approach. In the density- dependent simulation, the star for- 
mation rate slowly rises after periapse as gas interacts within 
the perturbed disks and is driven into central regions. The 
increasing central gas density eventually yields a broad peak 
in star formation rate at T c± —540 Myr. Due to orbital de- 
cay, the two galaxies return for a second passage and col- 
lide again at T ~ —260 Myr. This results in another sharp 
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Figure 5. Dynamical evolution and star formation in simulations of NGC 7252 in a progression of time steps equal to 155 Myr at 
T ~ —620, —465, —310 and —155 Myr (from left to right). Top Row: Gray scales show the gas distribution overlaid with contours as the 
distribution of old disk stellar populations. No star formation is included. Middle and Bottom Rows: In each image old stellar distribution 
is shown in halftone while star formation is shown as points in the density- dependent (middle) and the shock-induced simulations (bottom). 
Red points are stellar particles that have ages r < 10 Myr, green have r < 40 Myr, and blue have r < 155 Myr. 



peak in the shock-induced simulation, where fresh gas con- 
stantly returns to galaxies from the tails and settles in the 
disk plane. In the density- dependent simulation, following 
the pattern of the first passage, the star formation rate does 
not increase significantly until the gas is driven in towards 
the center. 

At T ~ —215 Myr the two galaxies merge. In the 
density- dependent simulation a sudden increase in gas den- 
sity in the nuclear regions results a strong and somewhat 
prolonged burst of star formation rate, which lasts to T c± 
— 155 Myr. In contrast, the shock-induced simulation shows 
a much more prompt but brief starburst. 

Figure [5] shows the dynamical evolution of gas, old disk 
stars and interaction-induced star formation, in a progres- 
sion of time steps equal to 155 Myr. The top row shows the 
distribution of gas and the old disk stellar populations with- 
out any star formation, while the middle and bottom rows 
contrast star formation in the density- dependent and shock- 
induced simulations. At each time in the density- dependent 



simulation, most ongoing star formation, indicated by the 
red points, occurs at the central regions of the galaxies. On 
the other hand, at the first passage in the shock-induced 
simulation, the starburst happens along the interface where 
the two disks collide, and the star formation extends to the 
outskirts of the gas disks. At T ~ —465 Myr some star for- 
mation products (blue and green points) populate the gas 
poor regions beyond the tip of the lower tail. These popula- 
tions formed from the impinging gas, which is converted into 
stars so promptly that it still retains much of its momentum 
and so populates the gas poor regions. Overall, the star for- 
mation has a wider spatial distribution in the shock-induced 
than in the density- dependent simulation. This larger spatial 
extent of star formation in shock-induced model is a corollary 
of the earlier onset of activity in such mod els, since the gas 
is more widely distributed at earlier times (Barnes 2004). 

Figure [6] shows matches to the present configuration 
of NGC 7252 for the two star formation prescriptions. The 
spatial extent of past star formation in the simulations is 
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Figure 6. Best match of the simulation to the present configuration of NGC 7252. Left: Gray scales show the gas distribution overlaid 
with contours as the distribution of old disk stellar populations. Middle and Right: In each image old stellar distribution is shown in 
halftone and star formation is shown as points in the density- dependent (middle) and the shock-induced model (right). Red points are 
stellar particles that have ages r < 155 Myr, green have 155-310 Myr and blue have 465-620 Myr. Same number of points are plotted 
in both simulations. 



Table 2. Median radii (in kpc) for the burst populations. 





Bi 


B 2 


B 3 


B 4 


all 


Density- dependent 




0.13 




0.07 


0.64 


Shock-induced 


4.32 


1.32 


1.29 


0.17 


2.11 



reflected in the distribution of star formation products at 
the present time. Here the stellar particles that formed dur- 
ing the merger are color coded according to their ages. The 
blue particles have ages of 465-620 Myr, corresponding to 
the products of the starburst occurred at the first passage 
(T ~ —620 Myr). The green particles have ages of 155— 
310 Myr, and the red particles have ages less than 155 Myr, 
reflecting regions of more recent star formation. Note that 
in both simulations the same total number of particles are 
plotted. In both cases, most of the on-going star formation 
is now confined to the central regions of the remnant body. 
However products of star formation are found throughout 
the tails and the remnant, especially in the shock-induced 
simulation. The prompt starburst occurred at the first pas- 
sage in the shock-induced simulation produces a wide distri- 
bution of star formation throughout the whole galaxy, while 
in the density- dependent simulation the first passage did not 
induce a prompt burst of extended star formation, and the 
past star formation products are far less distributed outside 
the central body. 

To quantify the spatial distribution of new stellar par- 
ticles, we define four distinct starbursts in the shock-induced 
simulation, as shown in Figure[4l 'Bi' (ages of 590-655 Myr), 
'B 2 ' (515-580 Myr), 'B 3 ' (245-285 Myr), and 'B 4 ' (170- 
225 Myr). These bursts are associated with the maxima 
in star formation rate. In the density- dependent simulation, 
'B2' (455-560 Myr) is denned as the broad star forma- 
tion feature occurred after the first passage, and 'B4' (160- 
230 Myr) is related to the star formation peak at merging. 
Table [2] lists the median radii in kpc for the burst popula- 
tions, measured at the present time. The last column "all" in 
the tables refers to the stars formed throughout the simula- 
tions from —930 Myr to present. The center of the galaxy is 
computed by locating the the minimum of the potential well 
of the old bulge particles. In the density- dependent simula- 
tion, most burst populations, and indeed most stars formed 
throughout the encounter are concentrated within the cen- 



tral kpc. On the other hand, the burst population formed 
at the first passage in the shock-induced simulation extends 
more than 4 kpc from the center. Overall, the radial extent 
of stars formed in the shock-induced simulation is nearly 4 
times larger than the extent of those formed in the density- 
dependent case. The total initial gas mass of each simulation 
is 9.4 x 10 9 Mq. At present, the density- dependent simula- 
tion produces, in total, 9.1 x 10 9 Mq of stellar particles and 
shock-induced simulation produces a comparable amount of 
8.7 x 10 9 M . 



4.2 Distribution of Star Clusters 

SS98 identified and measured ages of 7 proto-globular clus- 
ters and 1 young stellar association in the remnant. Their 
IDs and ages are compared with the simulated star forma- 
tion histories of NGC 7252 in Figure Among the globular 
clusters, six of them lie between projected distances of 3- 
11 kpc from the center of NGC 7252 and have ages of ~ 400- 
600 Myr, indicating that they formed early in the recent 
merger. Some of the cluster ages have a wide uncertainty; 
both of our simulations successfully reproduce the range of 
cluster ages, although clusters W6 and SI 05, with ages of 
580 and 600 Myr, are more consistent with the prompt star- 
burst at first passage in the shock-induced simulation. 

The spatial distribution of the observed clusters, how- 
ever, strongly discriminates between our models. To com- 
pare with the observed locations of these clusters (3-11 kpc 
from the center) , the bottom panel of Figure [3 shows the 
age distribution of stellar particles located within this an- 
nulus in our simulations. The density- dependent simulation, 
indicated as the grey line, produces a gradually declining 
distribution of ages, with almost all interaction-induced star 
formation within the very central regions. In contrast, the 
shock-induced simulation produces many more stellar par- 
ticles in this annulus and a sharp peak around the first 
passage ~ —620 Myr, showing that star formation occurs 
in more dispersed regions away from the center and that a 
high portion of stellar particles within this annulus formed 
during the starburst in the first passage. This result ap- 
pears consistent with the distribution of ages observed in 
SS98, suggesting that shocks can be an important trigger of 
for mation of the s e clus ters. 

iMiller et al.1 ([l997) found that the youngest stellar pop- 
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T (Myr) 

Figure 7. Comparison of age distribution of stellar populations. Top Panel: Global star formation history (in Mq yr _1 ) of NGC 7252 
shown from —900 Myr to present (T = 0). Color indications are the same as Figure 3] Cluster ages from SS98 are plotted as dots with 
their uncertainties. Cluster S101 has an age upper limit of 10 Myr and cluster S114 has an age of ~ 1 Gyr. Note that cluster S105 and 
S114 each has another possible age of ~ 200 and ~ 40 Myr. Bottom Panel: Histograms of number of stellar particles formed in the 
simulation, located within 3-11 kpc from the center, measured at present time. 



illations in NGC 7252, with ages of ~ 10 Myr, lie within 
2.0 kpc in the central disk. In our simulations, 93% of the 
stellar particles with ages ^ 10 Myr in density- dependent 
and 84% in the shock-induced models are formed within 
2.0 kpc from the center. A further inspection of the central 
region is shown in Figure [8] The top row shows the view of 
the central disk of our best-fit model from the sky and the 
bottom shows the view from the orbital plane. Despite the 
inconsistent orientation of our disk mentioned in Section 3, 
the size is consistent with t he observed scale of about 2.5 kpc 
in radius ([ S chweizerl 1 1 9 8 2| ) . The red points in Figure [8] rep- 
resent the stellar particles that have ages ^10 Myr; in both 
density- dependent and shock-induced models most of these 
young stellar particles are formed within the disk. 



5 PHOTOMETRIC PROPERTIES OF 
NGC 7252 

Figure [9] shows the simulated B and U — B images of 
density- dependent and shock-induced models which best fit 
the present configuration. The B band images of both mod- 
els show a dust lane across the central regions. There are 
ripples and loops around the remnant in both models, and 
star forming regions in the tails are seen as brighter clumps 
in the images. Since U — B traces young and massive stars, 
the U — B images indicate where the most violent star for- 
mation is. 

Table [3] lists the observed and modeled Mb and col- 
ors of NGC 7252; magnitudes and colors computed without 
extinction are included in the parentheses. The colors from 
the shock-induced model seem a closer match to the observed 















10" 







Figure 8. View of the central disk of our best-fit model from the 
sky (Top) and the orbital plane (Bottom). Left: Gray scales show 
the gas distribution. Middle and Right: In each image the gas disk 
is shown in halftone and star formation is shown as points in the 
density- dependent (middle) and the shock-induced model (right). 
Red points are stellar particles that have ages ^10 Myr, green 
have 10 < r < 20 Myr and blue have 20 < r < 50 Myr. Same 
stellar particles are shown in both views for each model, and the 
same total number of points are plotted in both simulations. 

values. Integrated colors of the density- dependent model ap- 
pear to be redder than those of the shock-induced model. 
This can be considered a direct result of the star formation 
prescription since star formation tends to concentrate in re- 
gions of high gas density in the density- dependent model, and 
thus gives a higher extinction and much redder colors than 
the shock-induced model. For example, in density- dependent 
model extinction causes Mb to be 1.51 mag fainter and the 
colors to be > 0.2 mag redder. However, the total Mb of 
both models are still fainter than the observed values. 

Despite the good match to the morphology of our best- 
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Figure 9. Comparison of simulated £> and U — B images of NGC 7252 using density- dependent (Top) and shock-induced (Bottom) star 
formation rule. In the B images, grayscale represents surface brightness ranging from \±b = 20 (dark brown) to 32 (white) mag arcsec -2 . 
In the U — B images grayscale ranges from U — B = —0.1 (red) to 0.6 (blue). Contours represent ub with levels of 26.0, 28.0 and 
30.0 mag arcsec -2 . 



fit model, the tails in our model seem to be too faint com- 
pared to the observations. Table |4] lists B magnitudes and 
B — R colors for the tails from the observations and from our 
best-fit model. HGvGS94 mentioned that the tails are 0.2- 
0.5 mag bluer in B — R than the outer regions of the main 
body, however they did not give specific values for the tails. 
Averaging the B — R colors of the outer regions of the main 
body (region A through C and H through M in HGvGS94) 
gives a value of 1.3, so the average B — R color of the tails 
is about 0.8-1.1. The values from our models are obtained 
by measuring our simulated images within the same areas of 
the tails from HGvGS94. In both models, we find that the 
tails are fainter by about 2 mag in Mb than the observed 
value and have redder B — R colors. 

Many factors could contribute to this result and it is 
beyond our scope to investigate each in detail. Neverthe- 
less it seems possible that the main reason is a lack of star 
formation in the tails, as reflected in the redder colors of 
our simulations. The lack of star formation could be due to 
insufficient material in the tails, or to an incomplete descrip- 
tion of star formation in the tails. A simple check of the gas 
in the tails helps to distinguish between these possibilities. 

In our best-fit simulations, there is about 4 x 10 8 Mq of 
gas in the NW tail and about 1.3 x 10 8 Mq in the E tail (Ta- 
ble 2}; in round numbers, our modeled tails have ~ 20% of 
the gas observed in the real tails. This deficiency very likely 
occurs because our gas particles were initially distributed 
like disk particles. If we were to use a more realistic gas dis- 
tribution for our progenitor galaxies, i.e. typical late-type 
gas disks with radii extending to about twice those of the 



stellar disks, a much larger gas mass in the tails could be 
achieved. 

However, even if the gas content of our tails could be 
increased by a factor of ~ 5 to match the observations of the 
NW tail, the star formation associated with this gas content 
is still probably insufficient. In the NW tail, about 9 x 10 4 
and 2 x 10 5 Mq of stars were born in the past 100 Myr in 
the density- dependent and the shock-induced model respec- 
tively. Thus we estimate the present star formation rate as 
0.002 Mq yr _1 in the NW tail, assuming a constant star 
formation rate for the past 100 Myr. According to equation 
JTJ), if we scale the gas content in the modeled tails to the 
observed amount of gas, the maximum star formation rate 
would become about 0.01-0.02 Mq yr _1 . This increase of 
star formation rate would bring the total B-band luminos- 
ity in the NW tail to about 0.5-1 x 10 9 Lq for each model, 
assu ming a Salpeter init ial mass function and solar metallic- 
ity ([Leitherer et al.|[l999h . The resulting Mb would become 
— 16.4 and —17.3 mag for the density- dependent and shock- 
induced model respectively, which are still fainter than the 
observed magnitude of —18.1. 

The above calculations show that even if we match the 
observed amount of gas in both tails, our models still would 
not reproduce the observed brightness of NGC 7252's tails. 
In fact, the tails in our models contain enough gas to account 
for the observed tail luminosities, if we consider a constant 
star formation rate in the tails since the first passage. For 
example, the NW tail has sufficient gas to sustain a star 
formation rate of ~ 0.72 Mq yr _1 for the past ~ 620 Myr. 
The stars formed would produce a B-band luminosity of 
~ 6 x 10 9 Lq, increasing the magnitude of the NW tail 
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Table 3. Photometric properties (in mag) of NGC 7252. 





M B 


U-B 


B-V 


V-R 


Observed 


-21.22 


0.17 


0.66 


0.74 


Density- dependent 


-20.08 (-21.59) 


0.24 (0.02) 


0.70 (0.41) 


0.55 (0.39) 


Shock-induced 


-20.68 (-21.54) 


0.16 (-0.03) 


0.60 (0.42) 


0.50 (0.39) 



to Mb — — 19 mag, or about one magnitude brighter than 
observed. 

These results strongly indicate that our tails appear 
fainter and redder than the observation due to an incom- 
plete description of star formation. Clearly, the star forma- 
tion rules could be improved. In low density environment 
such as tails, density- dependent rules may underestimate the 
star formation rate based on its dependence of gas density. 
On the other hand, shock-induced rule might not be use- 
ful since the tails are expanding kinematically. One way to 
solve this problem is to modify the adopted star formation 
rules. For example, another mode of star formation could be 
added to the present form, in which the star formation rate 
is proportional to the gas density to the order of unity, so 
that p* oc C* pg 1 u m + Co p g . This additional term would 
then represent a finite c hance (= Co At, where At is each 
simulation time-step; see lBarnesll2004h for gas to form stars 
regardless of the circumstances in the simulations. 



6 SUMMARY 

We describe new simulations of the galactic merger 
NGC 7252 with gas dynamics and star formation included. 
In our models NGC 7252 is formed by the merger of two 
similar gas-rich disk galaxies which fell together with an ini- 
tial pericentric separation of ~ 1.9 disk scale lengths about 
620 Myr ago. Starting with plausible models for the pre- 
encounter disks, our simulation produce fairly well matches 
with the observed morphology and kinematics. We empha- 
size that rather than a perfect reproduction, our goal is to 
present a plausible dynamical model of NGC 7252, which 
provides insight into the star formation history and matches 
the photometric properties; this nonetheless represents a 
new and significant step beyond pure dyn amical mode ling. 

Following the study of NGC 4676 by Barnes! (|2004l ) and 
IChien et al.1 (2007), we consider two star formation mech- 
anism in our sim ulation and use the a ges o f young glob- 
ular clusters from Sch weizer Seitzerl (11 998) to constrain 
the trigger of star formation. IBarne 

r<2QQ4D suggested that 
density- dependent star formation rule gives an incomplete 
description of large scale star formation in interacting galax- 
ies. In our simulations of NGC 7252, density- dependent star 
formation happens mostly at gas rich regions concentrated 
in galactic centers throughout the simulation, whereas more 
products of past star formation in shock-induced model are 
distributed at wider projected distances and in tails. The 
total mass born in bursts are on the same order in each 
model, yet the median projected radius of the burst popu- 
lations is 0.6 kpc in density- dependent model and 2. 1 kpc in 
shock- induced model. This supports the suggestion in I Barnes] 
( 2004) that shock-induced star formation may provide a bet- 



ter match to the extended distributions of young clusters in 
NGC 7252 and other merger remnants. 

When compared with observations, the star formation 
histories of both models successfully reproduce the range of 
ages spanned by the observed clusters. Those that have ages 
of ~ 570-600 Myr, in particular, are likely to form during the 
starburst at first passage, and this seems consistent with the 
prediction of the shock-induced simulation. Within the pro- 
jected annulus where the observed clusters lie, the age distri- 
bution of stellar particles formed in the shock-induced sim- 
ulation shows a peak around such age, while in the density- 
dependent simulation the age distribution is flat and declin- 
ing. Again, almost all interaction-induced star formation oc- 
curred in the very central regions for the density- dependent 
model, and the shock-induced model is more likely to explain 
the concentration of ages of clusters observed in SS98. 

We also describe a method to incorporate simple stellar 
population synthesis models into our dynamical simulations 
in order to compare the predicted photometric properties 
with the observed values. Beyond the attempt of matching 
the morphology as in many previous simulations of interact- 
ing galaxies, we also try to match the photometric proper- 
ties. While the overall magnitudes and colors are consistent 
with the observed values, the tails seem to be fainter and 
redder than observed. We argue that this lack of star forma- 
tion, reflected by the redder colors, is due to an incomplete 
description of star formation, rather than insufficient mate- 
rial in the tails. An additional mode of star formation which 
is independent of large-scale conditions could help to solve 
this problem. 

Finally we discuss some future directions for making 
our star- forming models more realistic, which could help 
improve the modeled photometric properties. First, as dis- 
cussed in Section 12.41 our models do not incorporate the 
gradual mass return from evolving stellar populations. In 
order to incorporate population synthesis models to pro- 
duce luminosities, we have to make adjustments which are 
not entirely straightforward. Moreover, recycled gas might 
help boost the star formation rate throughout the last stages 
of the merger, improving the luminosities and colors of our 
models of NGC 7252. 

Another concern regarding the colors of the merger rem- 
nant tails comes from the colors of the progenitor disks. Disk 
galaxies generally become moderately but systematically 
bluer outward from the bulge, and this effect is caused by 
the age and metallicity gradients in the stellar populations as 
well as internal dust extinc tion (e.g. lWalterbos &; Kennicuttl 
Il987l : iBell fe de Jondl2000h . In order to improve our model- 
ing of the progenitor galaxies, we could include radial gra- 
dients of age and metallicity in the initial stellar disks. This 
could influence the photometric properties of the remnant 
tails, since they mostly come from the outer parts of the 
progenitor disks. With regard to our models of NGC 7252, 
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Table 4. Properties of the tails. 







NW tail 






E tail 






M B 


B-R 


H I Mass 


M B 


B-R 


H I Mass 




(mag) 


(mag) 


(10 8 M ) 


(mag) 


(mag) 


(10 8 M ) 


Observed" 


-18.09 


0.8- 1.1 


20.8±0.5 


-16.84 


0.8 - 1.1 


11.0±0.4 


Density- dependent 


-15.96 


1.32 


4.5 


-14.60 


1.37 


1.4 


Shock-induced 


-16.12 


1.25 


4.2 


-14.87 


1.27 


1.3 



"Values adopted or calculated from Hibbar d et al 1 (Il994l ). 



it might also help to improve the luminosities and colors of 
the tails. However, unless star formation in the outer regions 
and tails can be sustained throughout the merger process, 
the outer populations will fade fast after the encounter, re- 
ducing the effect of age gradients on tail properties. 
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